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1 Introduction 

The use of near-infrared light in functional medi¬ 
cal imaging has gained increasing interest during re¬ 
cent years. Techniques using near-infrared light for 
imaging can be used to study, e.g., hemodynamic re¬ 
sponses on the cortex. Brain activation measurements 
can be performed using, e.g., MEG or fMRI, which, 
however, require heavy instrumentation. In medi¬ 
cal near-infrared imaging, the instrumentation can be 
made quite small in size and the costs can also be kept 
relatively low, which makes portable devices and the 
use of the technique for long-time monitoring pos¬ 
sible. Due to no significant interference with bio- 
magnetic measurements, the method can be combined 
with other modalities. Therefore, near-infrared imag¬ 
ing is an interesting alternative in performing or veri¬ 
fying brain activation measurements. 

In the following, we first describe the principal idea 
of medical near-infrared imaging and spectroscopy. 
We then outline the different measurement techniques 
and some aspects concerning the propagation of near- 
infrared light in tissue, especially considering the use 
of the method for imaging. 

The developments in modeling and instrumentation 
in our laboratory are then presented. A light transport 
model is described, and an approach to image recon¬ 
struction based on linearization is outlined. The basic 
design of the frequency domain imager and observa¬ 
tion about its performance are presented. 

2 The method 

2.1 Near-infrared spectroscopy and imaging 

The use of near-infrared light for spectroscopy and 
imaging is based on the different wavelength depen¬ 
dent absorption and scattering properties of tissues. 
The wavelengths used typically vary from 600 to 
1000 nm, a range of wavelengths often referred to as 
the optical window of tissue. At these wavelengths, 
many tissues are relatively transparent. At wave¬ 
lengths less than 600 nm, absorption by hemoglobin 
is very strong, whereas the strong absorption by water 


prevents using wavelengths beyond 1000 nm. 

The principal idea in near-infrared spectroscopy and 
imaging is to measure the properties of transmit¬ 
ted and/or reflected near-infrared light on the surface 
of the tissue at some distance apart from the light 
source. In spectroscopic applications, the concentra¬ 
tions of absorbing molecules such as oxy- and de¬ 
oxyhemoglobin and cytochrome oxidase can be es¬ 
timated. For imaging purposes, several light sources 
and detectors are used, and information about the in¬ 
ternal distribution of the absorbers and the scatterers 
can be retrieved. The aim can be, e.g., to localize 
changes in blood circulation and in the concentrations 
of absorbing molecules. 

2.2 Measurement techniques 

The methods in near-infrared imaging can be di¬ 
vided into three classes: continuous wave [1], time- 
resolved [2, 3] and frequency resolved [4, 5] tech¬ 
niques. In all of the techniques, the principle of the 
measurements is the same. Light from a light source, 
usually a laser, is delivered on the surface of the tissue 
using optical fibers. At the measurement site, on the 
surface of the tissue, typically a few centimeters apart 
from the source position, light is collected by another 
optical fiber and delivered to the detector. 

Continuous wave systems use continuous light and 
measure the intensity of the transmitted or reflected 
light. For imaging purposes, however, it is often pre¬ 
ferred to have more information on the propagation 
of light through tissue inhomogeneities. This can 
be achieved by measuring the temporal behaviour of 
transmitted light, either in the time or in the frequency 
domain. 

In the time domain techniques, a short laser pulse 
is used as a probe signal, and the distribution of the 
time-of-flights of the photons is measured. From this 
so-called temporal point spread function (TPSF), data 
such as the mean time of flight or some higher mo¬ 
ments can be retrieved. Alternatively, one can use a 
light source whose intensity is modulated with a high 
frequency, typically between 50 and 800 MHz, giving 
rise to a photon density wave propagating in the tis- 



sue. The amplitude and the phase shift of this wave 
relative to the incident light are measured. Theoreti¬ 
cally, by employing a range of different frequencies, 
the TPSF could be reconstructed. In practice, usu¬ 
ally a single frequency is employed, and the phase 
shift can be used to estimate the mean time of flight 
of the photons. The advantage of the frequency do¬ 
main method is that it is less expensive to implement. 

2.3 Imaging using near-infrared light 

In many biological tissues, the propagation of near- 
infrared light is dominated by scattering. Photons en¬ 
tering the tissue are generally multiply scattered be¬ 
fore being absorbed or reaching the detector. In a sin¬ 
gle measurement, a relatively large volume between 
the source and the detector is probed. Due to this type 
of diffusion of photons image reconstruction in near- 
infrared imaging is a challenging problem. 

Image reconstruction belongs to class of inverse prob¬ 
lems. In optical imaging, the inverse problem consists 
of finding the distribution of optical parameters in tis¬ 
sue, based on a set of measurements on the boundary. 
For solving the inverse problem, one first needs to 
solve the corresponding forward problem, i.e. to con¬ 
struct a model for light transport, that can be used to 
calculate the measurement set on the boundary, know¬ 
ing the distribution of optical parameters. 

3 Modeling and image reconstruc¬ 
tion 

3.1 Model for light transport 

Based on the particle interpretation, light can be ex¬ 
amined as a stream of energetic particles, photons. If 
the photon density is interpreted as proportional to the 
scalar field for energy radiance, various differential 
and integrodifferential equations can be established 
based on energy conservation. In the field of opti¬ 
cal imaging, the most widely applied equation is the 
radiative transfer equation (RTE) [6, 7], which is a 
balance equation describing the change of energy ra¬ 
diance in time due to changes in energy flow. 

For model construction, simplifications to the RTE 
are usually required. A common approach in appli¬ 
cations related to medical optical imaging is to use 
the diffusion equation (DE) [8], which can be derived 
from the RTE with a so called PI-approximation [6] 
and some assumptions valid for strongly scattering 
media [9]. 


If the measurements are performed in the frequency 
domain, the forward model can be constructed using 
the frequency domain diffusion equation, which is ob¬ 
tained from the time domain diffusion equation via 
Fourier transform: 

- — 4> - V ■ K,(r)W-A;$ a (r)4> = q Q . (1) 
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where = $(r,ca) is the complex fluence aris¬ 
ing from the intensity modulated part qj(r, u>) of the 
source term, oj the angular modulation frequency, 
k = 1/(3 (n' s + /!„)) the diffusion coefficient, /j, a the 
absorption coefficient and fl s the effective scattering 
coefficient. At the measurement site r m on the bound¬ 
ary, the complex exitance T, i.e. the outward directed 
flux, is calculated as: 

f(r m ,w) = -«(r m )n- V$(r m ,cu), (2) 

where n is the outward directed normal to the surface. 
The complex exitance is then used to obtain the mea¬ 
surable quantities, the amplitude 

A(v m M = |f (r m ,w).|, (3) 

and the phase angle 

(p(r m ,u)) = argr(r m ,cu), (4) 

of the photon density wave. 

As the boundary condition, so called Robin bound¬ 
ary condition [10], which can be derived from the as¬ 
sumption that the total inward directed energy current 
at the boundary is zero, is used: 

4(r,cu)+2«(r)h- V$(r,w) = 0 Vr £ <90. (5) 

The light source is modelled as an isotropic point 
source at a depth of 1 /(J s below the surface of the 
medium [10]. 

For iterative reconstruction algorithms, a relatively 
fast and flexible light transport model is desirable. For 
the numerical solution of the model described above, 
the finite element (EE) method was used. A two di¬ 
mensional FE model with linear basis functions was 
implemented. 

3.2 Solution to the inverse problem 

A widely used class of methods to address the inverse 
problem are techniques based on linearization and it¬ 
eration. In these approaches, the non-linear model U 
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4 Instrumentation 



Figure 1: Reconstruction of a perturbation in the ab¬ 
sorption coefficient p„ based on the logarithm of the 
modulation amplitude. On the background, p a = 
0.25 cm -1 and //'. = 20 cm -1 , and inside a sphere 
of a radius of 0.4 cm, centered at (0.7, —0.7), p a = 
2.0 cm~ l . 


between the optical parameters p and measurement 
quantities y is linearized at point po : 

y = U(p 0 ) + U'(po)(p-po) = yo + J8p (6) 

In Eq. 6, J, the derivative of the set of measurement 
quantities with respect to the set of optical parame¬ 
ters, is a large Jacobian matrix, which is needed for 
the reconstruction scheme. The FE model for the di¬ 
rect problem can be used to calculate J numerically. 
Using linearization and T ik honov regularization, the 
solution of the inverse problem can be expressed as a 
minimization problem: 

P = sol min {|| y - (y 0 + J5p )|| 2 + a\\L6p\\ 2 } (7) 

where is a a regularization parameter and L a regu¬ 
larization matrix incorporating prior knowledge. For 
L, e.g. a unit matrix can be used. From the mini¬ 
mization scheme, the gradient Sp can be solved and 
the solution presented as p = po + dp. In Fig. 1, 
a reconstruction of a perturbation on the absorption 
coefficient in a homogeneous background using cal¬ 
culated data is presented. The reconstruction is based 
on the logarithm of the modulation amplitude and a 
single step of the gradient. 


4.1 General 

The instrument under development in our laboratory 
is based on the frequency domain technique. The light 
source is a laser diode, which has a wavelength be¬ 
tween 750 nm and 850 nm. Optical powers between 
5 mW and 40 mW can be used with a high modula¬ 
tion efficiency. The light is intensity modulated using 
a modulation frequency in the RF region; 100 MHz 
and 450 MHz have been tested. In principle, with cur¬ 
rent laser diodes, the modulation frequency can be as 
high as 800 MHz. However, at higher frequencies, it 
is very difficult to retain the signal-to-noise ratio and 
dynamic range available at 100 MHz. 

The light is guided into an optical fiber and into a PC- 
controlled fiberoptic switch, which directs the light to 
one of 16 source fibers. Fight from the detection posi¬ 
tion is guided with a fiber bundle to a system consist¬ 
ing of a shutter and an aspherical lens, which focuses 
the light into a photomultiplier tube (PMT). The de¬ 
tection system is designed so that 1-4 bundles (one of 
which is selected active at a time) can be used with 
each detector. 

The optical fibers are attached to the tissue with a 
grid, which is heated and bent to match the shape of 
the tissue. In cases where there is no hair under the 
probes, a prism attachment can be used as an alterna¬ 
tive. 

The signal from the PMT is amplified and mixed with 
a second RF signal, which has a slightly different fre¬ 
quency from the modulation signal. The frequency 
difference is continuously compared with the inter¬ 
nal reference of the digital lock-in amplifier that per¬ 
forms the phase and amplitude measurement. The fre¬ 
quency of the secondary RF oscillator is continuously 
adjusted so that the two frequencies match. 

4.2 Instrument performance 

At 100 MHz, the detection limit of the instrument is 
approximately 30 nV rms / s/Wz, which corresponds to 
a light level of 70 nW rms at 800 nm. A typical signal 
amplitude with a 3 cm separation on the adult head 
with light hair is 0.3 mV. The signal to noise ratio of 
the amplitude in a phantom measurement at this sig¬ 
nal level is about 700/\/Hz, which results in a phase 
precision of 0.05°. In a physiological measurement, 
the noise is higher due to physiological changes in 
the tissue and the softness of tissue. The phase drift 
due to the electronics is approximately 0.1°/h for the 






References 



Figure 2: Block diagram of the instrument. 


first hour after startup, and 0.02°/h after that. The 
detector and the laser diode produce a larger drift, but 
temperature stabilization is expected to improve this. 
Accurate phase measurement requires that only the 
phase shift of the photon density wave affects the 
measured phase. In practice, a high anode current in 
the photomultiplier tube causes amplitude-dependent 
changes in the measured phase. This is a subject of 
ongoing research. 

5 Conclusions 

Initially, the instrument will be used for topograph¬ 
ical cortical activation studies, and in the future, for 
optical tomography. In three-dimensional imaging, 
a larger dynamic range of the measured intensity is 
necessary so that measurements at a wide range of 
source-to-detector distances can be performed. 

For imaging, more rigorous modeling and reconstruc¬ 
tion methods are being developed. The FE light trans¬ 
port model and the reconstruction algorithm based on 
linearization form a basis for further research. 
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